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Abstract 

This  report  summarizes  C3q>eriraents  and  analysis  concerned  with  sandwich 
and  cast  composite  solid  propellant  combustion.  Ihe  ingredients  used  in  the 
ejiperiments  are  ammonium  pershlorate  as  the  oxidizer,  hydroxjfl.  terminated 
polybutadiene  as  the  binder  and  four  catalysts:  Harshaw  catalyst  CU-0202, 

ferrocene  and  iron  blue.  Cinephotomacrography  and  the  fuze  wire  tech- 
nique are  used  for  combustion  visualization  and  burn  rate  determination. 
Scanning  electron  microscopy  is  used  for  quenched  sample  visualization. 

Areas  investigated  are  a)  analytical  and  experimental  determination  of 
synergistic  catalytic  effects  in  sandwich  and  propellant  combustion,  b)  the 
loading  of  ferrocene  into  the  binder  at  the  molecular  level  and  its  effect 
on  sandwich  combustion,  c)  differential  scanning  calorimetry  of  catalyst 
laden  binder  and  d)  analysis  of  sandwich  deflagration. 
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I.  Introduction 

TJiis  report  summarizes  efforts  of  the  second  year  on  contract  ONR  No. 
N0014-67“0159“0016-  Prior  results  were  concerned  with  the  sandwich  config- 
uration ahd  the  individual  effects  of  the  catalysts  Harshaw  cateilyst  CU0202, 
Fe„0_,  ferrocene  and  iron  blue.  Ihis  year's  effort  consisted  of  both  sand- 
wich  and  cast  propellant  investigations  to  determine  a)  possible  synergistic 
catalytic  effects  whereby  for  the  same  totad  loading  of  catalyst  the  effect 
of  two  catalysts  wo'-dd  be  greater  than  either  one  alone  and  b)  the  behavior 
of  a unique  type  of  ferrocene  loading  at  the  molecular  level  into  R-45 
prepolymer.  Differential  scanning  calorimeter  studies  were  also  performed 
on  catalyst  laden  hydroxyl  terminated  polybutadriene  (HIPB)  binder.  Analysis 
was  conducted  to  explain  synergistic  catalytic  effects  eind  to  improve  on  a 
previous  model  of  sandwich  deflagration. 

The  reason  for  investigating  synergistic  effects  is  that  they  were  ac- 
cidentally discovered  on  the  prior  year's  program.  Incorporation  of  this 
effect  would  allow  the  propellant  chemist  more  flexibility  in  propellant 
formulation.  The  investigation  of  loading  of  ferrocene  at  the  molecular 
level  was  motivated  by  prior  conclusions  that  physical  loading  of  catalyst 
into  the  binder  is  a poor  way  to  obtain  catalysis.  The  differential  scan- 
ning calorimeter  work  was  motivated  by  prior  results  that  catalysts  did 
not  appear  to  change  the  pyrolysis  mechanism  of  the  binder;  a confirmation 
was  desired  in  other  than  a deflagration  situation.  The  analysis  was  re- 
quired as  an  aid  in  reasoning  the  significance  of  the  experimental  results 
and  to  see  to  what  degree  a sandwich  analysis  may  be  used  for  a cast  can- 
posite  propellant  deflagration  analysis. 


II.  Determination  of  Possible  Synergistic  Effects 
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a.  IVo- Dimensional  Composite  Solid  Propellants 

iVo-dimensional  composite  solid  propellant  sandwiches  were  used  to  in- 
vestigate the  optimum  location  and  the  extent  of  the  catalytic  action  of 

(1) 

four  possible  burn  rate  modifiers.  Both  the  location  and  the  con^jounds 
were  varied  independently.  During  that  investigation  it  was  determined 
that  if  one  catalyst  was  added  to  the  oxidizer  prior  to  compaction  of  the 
disk  and  another  effective  catalyst  was  added  to  the  binder,  there  was  a 
net  increase  of  the  burn  rate  over  that  of  the  sum  of  the  independent 
actions  of  the  catalysts.  This  combined  effect  is  denoted  as  a synergistic 
effect  on  the  burn  rate.  The  current  study  has  systematically  investigated 
this  positive  synergism  for  all  combinations  of  the  same  four  con^iounds. 

The  two-dimensional  composite  solid  propellant  sandwiches  were  pre- 
pared by  the  method  outlined  in  Reference  ? from  disks  of  poly crystalline 
oxidizer  and  layers  of  binder.  The  disks  were  pressed  from  crystalline  am- 
monium perchlorate  (ap)  in  a mold  at  30,000  psi.  The  disks  were  assembled 
into  multi-layered  sandwiches  by  binding  them  together  with  layers  of  hy- 
droxyl terminated  polybutadiene  (HTPB) . ThvC  thickness  of  the  binder  layer 
was  maintained  at  I50  >im  during  the  em’ing  process  by  using  Teflon  spacers 
and  spring  loaded  sample  clamps. 

The  four  compounds  investigated,  weie  Harsb.aw  catalyst  CuD202,  (CC)  , 
ferric  oxide,  (lO),  iron  blue,  ( IB),  and  ferrocene,  (f)  . Hai's haw  catalyst 
Cu02C^2  is  a cormrcrcially  avaliable  catalyst.  It  is  a mixture  of  approxi- 
mately 82^;  cupric  oxide,  CuO,  and  VJio  chromic  oxide,  CigO^.  It  is  used 
in  accelerating  both  oxidation  and  reduction  processes.  It  is  commonly 
referred  r,o  as  copper  chromite.  Ferric  oxide,  exhibits,  as  do  most 


transition  element  compounds,  a high  degree  of  catalytic  activity.  It  is 
used  primarily  for  oxidation  processes.  Iron  blue  is  commericaH^  used  as 
a pigment.  It  is  a complex  ammonium  iron  hexacyanoferrato  with  the  chemical 
formula,  Fe(KHj^)Fe(CW)g  . It  has  a cubic  crystalline  structure.  Ihe  iron 
is  present  as  both  ferric  and  ferrous  ions  in  the  lattice.  Ferrocene  is  an 
orgeuiometallic  compound  with  the  chemical  name,  biscydopentadienyl  iron, 
(C2H^)gFe.  It  is  a yellow  crystalline  solid  with  relatively  high  thermal 
stability  for  an  crganometallic  compound. 

Ihe  samples  were  prepared  with  constant  volumetric  loading  of  the 
compounds  based  on  an  addition  of  2v4  to  the  oxidizer.  Oils  amount  of 
catalyst  was  added  to  the  AP  prior  to  grij  '.ng  and  pressing  and  H.37w^  of 
catalyst  was  added  bo  the  HTPB. 

Ihe  cured  oxidizer-fuel  sandwiches  were  cut  into  8 mm  by  4 ram  samples. 
These  prepared  samples  were  mounted  in  the  pressurized  combustion  appara- 
tus  of  Jones.  Motion  pictures  of  the  burning  sample  were  obtained  at 
a rate  of  I6OO  frajnes  per  second  at  a latentimage  magnification  of  1:1  and 
2:1.  These  motion  pictures  were  used  to  obtain  an  accurate  value  of  the 
sample  burning  rate.  Ihe  sandwich  vertical  burn  rate  and  the  burn  rate 
normal  to  the  oxidizer  surface,  as  defined  in  Reference  (l) , were  obtained 
for  all  combinations  of  the  four  catalysts  along  with  suitable  base  line 
cases  at  the  same  volumetric  loading  in  the  binder  and  oxidizer  at  a com- 
bustion pressure  of  600  psia. 

Ihe  tabulated  results  for  the  determination  of  possible  synergistic 
effects  are  shown  in  Table  1.  The  base  line  comparison  cases  for  a single 
catalyst  present  were  for  of  catalyst  added  to  the  AP  and  ’:.37w^  of 
catalyst  added  to  the  HTPB.  This  gave  a uniform  volumetric  loading  of 


TABLE  I.  BURN  RATES  FOR  DETERMINATION  OF  SYNERGISTIC 
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catalyst  across  the  binder-oxidizer  interface  of  the  two-dimensional 
sandwiches.  The  total  weight  of  the  catalyst  was  maintained  constant  for 
the  conl)inations  of  catalysts.  One  weight  percent  of  one  catalyst  along 
with  Iw^  of  another  catalyst  was  added  to  the  AP  and  2.2w^  of  each  catalyst 
was  added  to  the  HTPB.  The  results  obtained  for  the  base  line  cases  of 
single  catalysts  added  to  the  sandwiches  are  in  agreement  with  previous 
data.^^^  It  must  be  noted  that  the  comparison  cannot  bo  exact  since  in 
the  previous  investigation the  catalyst  loading  was  discontinuous  across 
the  binder-oxidizer  interface,  but  in  previous  e:q>eriments  catalyst  in  llie 
binder  was  shown  to  be  ineffective. 

Ferric  oxide,  iron  blue  and  ferrocene,  when  used  alone,  were  burn  rate 
inhibitors  for  AP  at  600  psia.^^^  Ihe  burn  rate  normal  to  the  oxidizer 
surface  was  less  than  that  of  the  uncatalyzed  AP.  This  is  also  true  in 
Table  1.  When  these  proven  irihibitors  of  the  AP  deflagration  process  were 
used  in  coni)inations  with  the  same  total  mass  of  catalyst  present,  the  burn 
rate  normal  to  the  oxidizer  surface  was  increased  over  that  of  the  un- 
catalyzed AP.  This  was  noted  for  the  10- IB,  lO-F  and  IB-F  systems .Harshaw 
catalyst  Cu0202  increased  the  AP  deflagration  rate  when  used  alone  at  600 
psia.  All  combinations  of  Harshaw  catalyst  Cu0202  with  the  other  three 
compounds  exhibited  a synergistic  effect  on  the  AP  deflagration  rate  even 
thoi^gh  three  of  them  acted  as  inhibitors  when  used  alone. 

The  maximum  saridwich  burn  rate  (.5^  inches/sec)  has  been  obtained  witJa 
the  combination  of  PEarshaw  catalyst  Cu0202  and  ferric  oxide.  The  largest 
synergistic  effect  (burn  rate  change  as  compared  with  either  single  cata- 
lyst alone)  has  been  seen  for  a ferric  oxide  and  iron  blue  mixture  with  a 
burn  rate  of  .tl  inches/sec.  All  combinations  of  catalysts  did  give  some 
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degree  of  positive  synergism  of  the  sandwich  vertical  burn  rate.  For  the 
combinations  with  Harshaw  catalyst  Cu0202  only  ferric  oxide  was  able  to 
clearly  increase  the  burn  rate  above  that  of  Harshaw  catalyst  CuD202  alone. 
All  combinations  of  ferrocene  in  this  coni^iguration  yielded  relatively  low- 
er santh.’ich  vertical  burn  rates. 

Considering  these  results,  the  most  promising  two  combinations  of 
catalysts,  Harshaw  catalyst  Cu0202  - ferric  oxide  and  iron  blue  - ferric 
oxide  were  chosen  for  testing  over  the  pressure  range  of  300  to  2000  psia. 

Ihe  results  for  the  Harshaw  catalyst  Cu0202-ferric  oxide  system  tests 
are  shown  in  Figure  1.  The  maximum  effect  is  at  1500  psia.  Both  catalysts 
were  equally  effective  in  catalyzing  the  burn  rate  when  used  alone.  When 
half  of  one  catalyst  was  replaced  by  the  same  amount  of  the  second  catalyst 
the  burn  rate  doubled.  The  sandwich  vertical  burn  rate  and  the  burn  rate 

>:  normal  to  the  oxidizer  surface  both  exhibited  positive  synergism  over  the 

i 

entire  pressure  range. 

For  the  iron  blue  - ferric  oxide  samples  the  results  were  not  as  clear. 
They  are  shown  in  Figure  2.  The  maximum  positive  effect  occurs  at  600  psia  . 
It  is  not  conclusive  at  1000  and  I500  psia  that  a positive  synergistic  ef- 

i 

feet  has  been  obtained,  since  the  combined  catalysts  produce  a burn  rate 

j that  lies  between  the  burn  rates  obtained  for  a single  catalyst  present. 

\ A negative  effect  has  not  been  produced,  i.e.  the  burn  rates  have  not  been 

♦ 

i 

! inhibited.  Within  the  experimental  accui-acy  of  the  test  the  burn  rate  nor- 

mal to  the  oxidizer  surface  for  the  combined  catalysts  loading  has  followed 
i the  burn  rate  for  ferric  oxide. 

As  in  previous  tests  the  separation  between  the  dashed  and  solid  curves 


is  representative  of  the  amount  of  catalytic  activity  talcing  place  in  the 


C«i0202 
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PRESSURE  (PSIA> 

FIGURE  I.  BURN  RATE  RATIO  FOR  HARSHAW  CATALYST  Cu  0202  AND  FERRIC  OXIDE. 


1%  IRON  BLUE  a 1%  FERRIC  OXIDE 


PRESSURE  (PSIA) 

FIGURE  2.  BURN  RATE  RATIO  FOR  IRON  BLUE  AND  FERRIC  OXIDE 
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binder-oxidizer  reactions.  The  tests  that  showed  large  positive  syner- 
gistic effects,  Harshaw  catalyst  Cu0202  and  ferric  oxide  at  1500  and  2000 
psia,  and  iron  blue  and  ferric  oxide  at  600  and  2000  psia,  also  showed  sub- 
stantial increases  in  this  separation  of  the  two  burn  rates. 

b.  Cast  Composite  Propellants 

Maintaining  the  same  volumetric  loading  of  catalysts  as  for  the  two- 
dimensional  sandwiches,  three-dimensional  cast  composite  solid  propellant 
strands  were  prepared  to  see  if  this  synergistic  effect  would  be  carried 
over  to  the  real  cast  propellant.  These  strands  were  prepared  from  a com- 
mon lot  of  uncured  composite  propellant  with  an  solids  to  binder 

loading.  The  oxidizer  was  chosen  to  have  a bimodal  ammonium  perchlorate 
particle  distribution  of  30^  UO  and  70?&,  iSC  pjn.  The  binder  was  hydroxyl 
terminated  polybutadiene.  Burning  rate  data  are  available  for  this  un- 
catalyzed propellant . It  was  similar  to  propellant  ^8  in  the  Princeton 
University  test  series.  The  catalyst  was  added  to  the  unciu-ed  propellant 
at  a weight  percent  of  2.41.  This  would  give  the  same  equivalent  volumetric 
loading  as  in  the  sandwiches. 

The  strands  were  cast  in  teflon  molds  and  all  samples  were  cured  for 
the  same  length  of  time.  The  strand  size  was  .25  inch  by  .2'^  inch  by  1.7 
inches  long.  The  burn  rates  were  obtained  by  a fuse  wire  technique  ' 
using  seven,  ^ ampere  (~  .010  inch  diameter) , wires  inserted  in  predrilled 
holes  which  were  spaced  .2  inches  apart.  Each  fuse  wire  was  connected  in 
series  to  a 10  KH  resistor.  These  seven  resistors  were  connected  in  par- 
allel with  a 57  Kp  resistor  across  an  18  volt  battery  power  supply.  As 
each  fuse  wire  burns  through  the  equivalent  resistence  increases  and  the 
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voltage  is  recorded  on  an  oscillograph. 

The  strand  with  fuse  wires  is  shown  in  Figure  3-  It  is  mounted  on 
the  has  of  the  combustion  bomb.  A sufficient  nitrogen  flow  is  passed 
over  the  sample  to  prevent  recirculation  of  the  hot  gases  that  can  cause 
premature  melting  of  the  fuse  wires.  This  nitrogen  flow  has  varied 
from  .7  to  1.4  ft/sec.  ihe  most  stable  resxilts  were  obtained  at  the 
higher  flow  rates . Sample  oscillograph  traces  are  shown  in  Figure  4 for 
an  uncatalyzed  strand  and  one  with  ferric  oxide  added  to  the  strsind. 

The  burn  rates  for  the  addition  of  the  two  most  promising  combi- 
nations of  catalysts  as  determined  from  the  two-dimensional  composite 
solid  propellant  test  are  tabulated  in  Table  2.  These  strands  of  pro- 
pellant were  burned  at  6OO  psia.  For  the  combination  of  ferric  oxide 
and  iron  blue,  the  burn  rate  is  between  the  burn  rates  of  the  samples 
with  only  one  catetlyst  present.  The  ferric  oxide  - Harshaw  catalyst 
Cu0202  system  showed  a reduction  of  burn  rate  below  that  of  either  of 
the  two  catalysts  used  alone. 

The  burn  rates  do  not  compare  f?’«forably  with  the  res\ilts  obtained 
in  Table  1 for  the  two-dimensional  sandwiches.  The  burn  rate  for  the 
uncatalyzed  strands  agree  with  the  Princeton  data.  The  stvuctural  strength 
and  elasticity  of  the  samples  containing  both  Harshaw  catalyst  Cu0202 
and  iron  blue  were  considerably  different  from  the  ferric  oxide  and  pure 
samples.  The  ferric  oxide  and  pure  sainples  were  softer  and  more  easily 
removed  from  the  muH  v^iilo  the  samples  containing  iron  blue  and  Harshaw- 
catalyst  CuD202  were  brittle  and  exhibited  several  tests  with  very  rapid 
burning.  These  tests  were  not,  considered  in  the  preparation  of  Table  2. 


FIGURE  3. 


A STRAND  WITH  SEVEN  FUSE  WIRES  MOUNTED  ON  THE  COMBUSTION 


BOMB  BASE 


TABLE  H.  BURN  RATES  FOR  CAST  COMPOSITE  PROPELLANT  STRANDS. 
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These  cast  propellant  samples  were  initially  prepared  to  contain  on 
average  the  same  volumetric  loading  of  catalyst  as  tlie  sandwiches.  But 
for  the  cast  propellants  all  of  the  catalyst  is  suspended  in  the  binder 
matrix,  which  constitutes  I7v%  of  the  sample.  Since  2.4lw^  of  catalyst 
was  added  to  the  cast  propelleint  samples,  the  binder  ms-trix  contains 
14.2w^  of  the  catalyst.  The  binder  matrix  for  the  sandwiches  contained 
4.4^  of  catalyst.  Uiis  increase  of  catalyst  in  the  binder  effected  the 
cure  and  the  strength  of  the  binder.  This  high  catalyst  loading  in  the 
binder  may  be  above  the  amount  required  for  maximum  b;irn  rate  augmen- 
tation. It  was  determined  eai'lier^^^  that  the  addition  of  the  catalyst 
directly  into  the  polycrystalline  structure  was  more  effective  than  addi- 
tion to  tne  binder.  If  the  catalyst  could  be  distributed  in  the  oxi- 
dizer crystalline  structure  to  give  the  same  volumetric  loading  as  in 
the  sandwiches  a more  realistic  comparison  could  be  made.  This  would 
require  a developm.ent  of  a n_w  ajnmonium  perchlorate  production  technique. 

It  may  be  possible  to  establish  the  existence  of  a synergistic  effect 
in  cast  propellants  if  the  total  catalytic  loading  in  the  binder  is  re- 
duced to  a more  compardble  level  with  the  two-dimensional  sandwiches. 

Tills  approach  will  be  followed  in  future  efforts. 

III.  Indian  Head  Ferrocene 

It  had  been  a prior  conclusion^^^ ’ that  crystalline  catalysts 
physically  loaded  into  the  binder  were  ineffective  in  augmenting  sand- 
wich deflagration  rates.  The  apparent  reason  is  that  with  the  catalyst, 
in  the  binder  the  diffusion  of  crystalline  matter,  once  it  is  released 
into  the  gas  phase,  toward  the  hot  AI’  is  too  slow  to  augment  the  binder-oxidizer 


15 


reactions  fast  enough  to  increase  heat  i'cedhack  to  the  sandwich  surface. 
Furthermore,  the  dlspers'on  of  catalytic  material  must  be  greater  for  a 
substance  dispersed  at  the  molecular  level  rather  than  at  a macroscopic 
crystalline  level.  Following  this  hypothesis,  it  ’nDuld  be  desirable  to  in- 
troduce the  catalys . in  molecular  rather  than  crystalline  form,  once  the 
binder  had  pyrolysed  to  the  gas  phase.  The  diffusion  of  gas  phase  molecules 
should  be  faster  than  crystalline  particles. 

It  was  learned  that  the  Naval  Ordnance  Station  at  Indian  Head, 

Maryland  was  working  in  the  area  of  chemically  loading  ferrocene  into  R-h5 
polymer. Indian  Head  supplied  some  of  the  polymer  with  Fe  by  weight 
in  the  pol3cner.  Due  to  an  advanced  state  of  cure  of  the  recr-ived  material 
work  with  this  substence  was  limited,  but  one  sandwich  was  prepared  with 
this  material. 

Figure  5 shows  scanning  electron  microscope  pictures  of  quenched  sam- 
ples burned  at  1000  psia.  One  sainple  shown  is  from  prior  work  with  the 
crystalline  ferrocene  in  the  binder  alone  and  the  other  sample  is  with  the 
Indian  Heeid  ferrocene.  The  striking  effect  of  the  chemically  loaded  ferro- 
cene is  the  depression  of  the  AP  in  tnc  vicinity  of  the  oxidizer -binder 
interface.  This  is  interpreted,  as  shown  by  many  past  results,  as  increased 
chemical  rates  between  the  binder  and  oxidizer  for  the  Indian  Head  binder  as 
compared  with  the  binder  with  cr;/stallinc  catalyst.  Because  the  AP  is  not 
flat  (horizontal)  the  vertical  burn  rate  of  the  sandwich  had  to  be  greater 
with  the  Indian  Head  binder.  This  result  is  considered  to  confirm  the 


hypothesis  of  molecular  diff’asion  and  greater  dispersion  as  explained  above. 
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IV,  Differential  Scanning  Calorimetry 


Since  prior  sandwich  results  led  to  tlie  conclusion  that  catalytic 
effectiveness  is  poor  when  the  catplyst  is  loaded  into  the  binder,  and 
there  has  been  no  evidence  that  catalysts  increase  the  binder  j^Tolysis 
rate,  it  was  desired  to  check  these  conclusions  by  differential  scanning 
colorimetry  applied  to  catalyst -laden  binder.  The  apparatus  used  was  a 
Perkin-Elmer  Differential  Scanning  Calorimeter,  Model  DSB-lB,  and  the 
procedures  used  were  identical  to  those  of  Ref.  6.  Ihe  binder  used  was 
HIRE  and  two  catalysts,  iron  blue  ind  ferrocene,  were  investigated.  The 
nominal  weight  of  each  sample  was  5 mg.  and  the  amount  of  catalyst  load- 
ing was  consistent  with  that  used  in  the  sandwich  deflagration  tests. 


Tiie  DSC  unit  measures  the  difference  in  power  required  to  heat  two 


samples,  one  inert  and  one  containing  the  decomposing  sample,  through 
a programmed  temperature  rise.  Sample  graphs,  shown  on  Figure  6 are  two 
UTP3  runs  at  a scan  rate  of  5 K/min.  ihis  figure  shows  the  reproduc- 
ibility attai^ble  between  runs.  2ie  heat  of  gasification  is 

. 

q = j Q dt  = , Q (l/dT/dt)  dt  where  T„  is  the  teraperatiu'e  at  iriiich  gasi- 

C T 

fication  for  puie  HTPB  are  shown  in  Table  3*  Ih  Table  3 there  is  a high 
degree  of  data  scatter  but  an  vinmis takeable  trend  of  heat  of  gasifi- 
cation with  scan  rate  depe'ident.  A plot  of  q v.s.  l/scan  rate  is  given 
in  Figure  7.  A line  through  the  most  consistent  set  of  data  (HIPB-IB) 
shows  an  extrapolation  to  infinite  scan  rate  of  q = 26^?  cal/g.  TJiis 
should  be  compared  with  433  cal/g  obtained  in  Ref.  7>  vrtiich  was  obtained 
at  heating  rates  comparable  to  those  found  during  propellant  deflagration. 


FIGURE  6.  DIFFERENTIAL  HEAT  INPUT  AS  A FUNCTION  OF  TEMPERATURE. 


Table  3 

Heat  of  Gasification  for  HTPB  and  HTPB  - catalyst 


Rate  (K/min) 

Binder 

No.  of  Tests 

q (cal/g) 

20 

HTPB 

1 

328 

10 

HTPB 

3 

518  * 

5 

HTPB 

2 

+ 129 
735  - 129 

20 

HTPB-IB 

1 

345 

10 

HTPB-IB 

5 

^15 : % 

5 

HTPB-IB 

1 

584 

20 

HTPB-F 

1 

285 

which  are  orders  of  magnitude  higher  than  may  be  found  with  the  DSC  unit. 

The  primary  conclusion,  however,  is  that  within  the  data  scatter  there  ap- 
pears no  definite  effect  of  IB  or  F on  the  heat  of  decomposition  of  HTPB. 

Figures  8 and  9 show  the  actual  traces  for  pure  HTTO  and  HTPB-IB  at 
scan  rates  of  10  K/min.  The  catalyst  has  very  little  effect  upon  the  tem- 
peratxire  at  the  onset  of  decomposition  activity.  There  does  appear,  however, 
to  be  a minor  effect  of  IB  in  smoothing  out  the  roughness  near  the  peih  of 
the  decomposition  activity.  Nevertheless,  the  overall  temperature  width 
of  activity  and  the  heat  of  gasification  is  affected  in  only  a minor  way. 

The  same  conclusion  »”ay  be  derived  from  the  HTPB-F  run. 


V.  Ansilytical  Efforts 

IVo  major  analytical  efforts  have  been  conducted  during  the  current 
year.  One  task  concerned  an  explanation  of  synergistic  effects  and  the 

/g\ 

other  concerned  an  extension  of  a previous  sandwich  deflagration  analysis . ' ^ 


During  the  course  of  extension  of  the  sandwich  analysis,  two  major  items 
were  discovered:  a)  the  solution  procedure  used  in  the  previous  effort  was 
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FIGURE  8.  DIFFERENTIAL  HEAT  INPUT  AS  A FUNCTION  OF 
TEMPERATURE. 
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FIGURE  9.  DIFFERENTIAL  HEAT  INPUT  AS  A FUNCTION  OF  TEMPERATURE. 
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in  error  eind  b)  an  inproved  solution  could  be  foxmd  for  the  solid  phase 
heat  transfer  problem.  As  a consequence  it  was  found  that  the  sandwich 
burn  rate,  which  could  not  be  found  with  the  old  procedure,  appears  as  ein 
eigenvalve  to  the  problem.  Ihe  restricted  problem  of  an  inert,  dry  binder 
adjacent  to  deflagration  AP  was  therefore  solved  by  the  new  procedure.  The 
details  are  located  in  the  Appendix. 

5yiiergistic  effects  were  first  treated  by  insertion  of  two  catalytic 
reactions  in  the  chemical  kinetics  scheme  of  Ref.  9 for  deflagrating  AP. 

It  is  believed  that  the  details  should  not  be  presented  here  because  it 
may  be  construed  that  these  catalytic  reactions  are  being  proposed  as  actual 
catalytic  mechanisms;  vdiereas,  they  were  only  introduced  as  a plausibility 
argument  for  synergistic  effect.  The  plausibility  argument  which  is  pre- 
sented here  encompasses  all  of  the  results  that  were  obtained  from  the  AP 
deflagration  model  without  confusing  the  issue  with  AP  details. 

To  explain  the  results  of  the  AP  catalysis  model  on  a simple  basis  the 
effect  of  two  catalysts  on  the  reaction  rate,  R,  may  be  modeled  as  follows: 


or 


"a  = "o 


(1) 

(2) 


^rtiere  the  f's  and  g’s  as  functions  of  the  catalyst  mass  fraction,  y^,  are 
shown  in  Figure  10.  Equations  (1)  and  (2)  merely  represent  two  functional 
assumptions  >riiich  yield  expected  behavior  with  respect  to  each  catalyst. 
There  is  a general  rise  in  rate  with  catalyst  loading  to  a maximum,  and  a 
fall  off  in  rate  to  zero  must  occur  as  the  catalyst,  \diich  is  in  inert, 
saturates  the  propeilemt.  Since  the  burn  rate  is  roughly  proportional  to 


26 


the  square  root  of  the  reaction  rate^^^  Eqs.  (1)  and  (2)  yield  the  two 
options  for  the  hui*n  rate: 


(3) 


where  K is  some  const8int,As  long  as  the  shape  of  the  rate  curves  in  Figures 
10a  and  10b  is  as  shown,  the  following  conclusions  are  valid: 

a)  r(Y  , ) > r(Y  , O)  or  r(0,  Y ) for  Y s Y 

'^l  ^2  ^1  ^2  ^ '^max 


and  b) 


r 

max 


max  ‘■max 


In  this  case  the  addition  of  a second  catalyst  will  always  augment  the  rate 
above  what  may  be  accomplished  with  a single  catalyst.  Moreover,  for  a 
fixed  total  mass  of  catalyst. 


+ Y = 
^2 


c, 


and  for  equally  effective  catalysts  so  thot 

fl  - fg  or  . gg 

r(Y  , c - Y ) > r(c,  0)  or  r(0,  c)  , 

Cl  c^ 


as  long  as 

0 < - Y 

®1  ^2 
max  max 

therefore,  with  the  appropriate  restrictions,  this  simple  demonstration 
shows  that  the  presence  of  one  catalyst  augments  the  effect  of  the  presence 
of  another  catalyst.  Ihis  is  a synergistic  effect  upon  rate.  Of  course, 


this  is  only  a plausibility  argument  and  no  claim  is  being  made  that  this 
is  the  actual  mechanism  or  that  Eqs.  (l)  and  (2)  are  generally  valid. 
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VI.  Conclusions 

1.  All  possible  combinatiuns  of  the  four  catalysts  investigated  ex- 
hibited a positive  synergistic  effect  on  the  sandvich  vertical  bur-n  rate 
when  tested  at  600  psia  and  a constant  volumetric  loading  based,  on  a total 
of  2w^  addition  of  catalyst  to  the  oxidizer. 

2.  Ihe  maximum  sandwich  burn  rate  at  600  psia  was  obtained  for  the 
combination  of  Harshaw  cateilyst  CuD202  and  ferric  oxide.  Vhen  these  samples 
were  investigated  over  the  pressure  range  of  300  to  2000  psia,  the  maximum 
synergistic  effect  occurs  at  1500  psia.  There  was  a positive  synergistic 
effect  over  the  entire  pressure  range. 

3.  The  maximum  synergistic  effect  with  the  two-dimensional  sandwiches 
at  600  psia  was  obtained  for  the  iron  blue  and  ferric  oxide  combination. 

This  combination  was  not  as  effective  as  the  Harshaw  catalyst  CuC'202  and 
ferric  oxide  over  the  entire  pressure  range. 

4.  The  cast  composite  propellant  strands  did  not  exhibit  a positive 
synergistic  effect  when  tested  at  600  psia. 

5.  It  is  suspected  that  the  heavy  catalyst  loading  is  responsible 
for  the  absence  of  synergistic  effects  in  the  cast  propellants.  A lower 
overall  loading  should  be  investigated. 

6.  From  testing  of  ferrocene  loaded  into  R-45  polymer  at  the  molecular 


level  it  is  concluded  that  molecular  loading  of  a catalyst  into  a propel- 
lant is  superior  to  physical  loading  of  a crystalline  material,  as  sus- 
pected from  previous  sandwich  testing. 
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7.  Differential  scanning  calorimetry  of  iron  "blue  and  ferrocene -loaded 
HTPB  confi  rms  earlier  conclusions  that  these  catalysis  do  not  modify  the 
pyrolysis  mechanism  of  the  hinder  to  an  appreciable  extent. 

8.  For  dry  inert  sandwiches  analysis  indicates  that  in  the  pressure 
range  of  20  to  100  atm  the  sandwich  should  burn  with  a nearly  flat  surface 
and  very  close  to  the  A?  deflagration  rate,  as  has  always  been  observed  for 
uncatalysed  sandwiches  with  even  wet  binders. 

9.  A sharp  slope  discontinuily  is  predicted  at  the  binder-oxidizer 
interface  for  dry  binders,  as  has  been  seen  experimentally  for  catalys'^d 
sandwiches . 

10.  For  typical  particle  sizes  of  the  order  of  20  lun  and  larger  the 
sandwich  analysis  can  form  the  basis  for  a composite  propellant  deflagration 
theory,  if  the  binder  is  dry. 
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ABSmCT 

i 

A theoretical  solution  has  been  obtained  for  the  shape  of  a deflagrating 
aaajonium  perchlorate  surface  ^en  it  is  adjacent  to  an  inert,  pyrolysing,  dry 
binder,  Die  eigenvalue,  the  regression  rate,  is  shown  to  be  independent  of 
binder  type  and  very  close  to  the  burn  rate  of  pure  AP,  as  has  been  ejqjeri- 
aentally  observed.  A slope  discontinuity  of  the  surface  should  exist  at  the 
binder-oxidizer  interface  and  typical  binders  should  incline  very  near  to  go° 
to  the  nearly  horizontal.  AP  surface,  at  the  junction  of  the  two.  !Ihe  transi- 
tion from  AP  to  binder  should  take  place  on  a distance  scale  of  the  ordei  of 
microns  and  all  results  are  quite  insensitive  to  pressure  level. 


This  work  was  sponsored  by  the  Office  of  Naval  Research,  Power  Branch, 
under  Contract  No.  KOOOl4-67-Ol59"^^l6-  Substantial  aid  in  the  calculations 
was  given  by  Mr.  Narendra  Kumar. 

Index  categories:  Coiribustion  in  Heterogereous  Media,  Solid  and  l^brid 

Rocket  Engines. 

Professor,  School  of  Aerospace  Engineering,  Member  AIAA. 
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g + Sqj^Yj, 

pyroJysis  law  constant  of  vapor  pressure  constant 
dimensionless  pyrolysis  law  constant  or  vapor  pressure  constant 
solid  phase  specific  heat 

specific  heat  at  constant  pressure  for  gas  phase 

flaae  standoff  distance 
deviation  of  c from  planar  AP  case 
activation  energy 

T/To 

deviation  of  temperature  from  planar  AP  case 
solid  phase  temperature  perturbation  constants 

preejqponential  fantor  in  reaction  rate  law 
dimensionless  preexponential  factor  in  reaction  rate  law 
constants  in  eigensolution 

coordinate  normeU.  to  the  solid-gaua  interface  directed  toward, 
the  gas  phase 

pressure 

reaction  integral  defined  by  Equations  (18) 

exothermic  reaction  heat  in  gas  phase,  endothermic  in  solid  phase 
xmiversal  geus  constant 
Reynolds  number 

burn  rate  auid  regression  rate  normal  to  surface,  respectively 
ten5)erature 

gas  velocity  in  y direction 
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Subscripts 

B 

f 

F 

g 

0 
s 
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production  rate  of 
coordinates 

perturbed  y position  of  solid  surface 
mass  fraction  of  species  k 

deviation  of  mass  fraction  from  planar  AP  case 

1 + (dy^/dx)^ 

thermal  diffusivity,  X/pc 

inclination  angle  of  binder 

dimensionless  activation  energy,  E/rt 

c /c 
s'  p 

thermal  conductivity 

c X /c  X 
p s s g 

density 

inclination  angle  of  solid  AP  surface 
binder 

flame  temperature 

NH3 

gas  phase 
cold  solid 

solid  phtse  or  surface 

quantity  evaU-uated  at  flame  standofi  position 


Superscripts 


quantity  cvailuated  for  the  one-dimensional  AP  deflagration 


/ 
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A 


ordinary  derivative  with  respect  to  x 

dimensional  quantity 

solid  phase  perturbation  function 
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A substantial  number  of  experimental  investigations  of  composite  solid 
propellant  ingredient  behavior  have  dealt  with  the  sandwich  configuration 
(1-8) . By  a sandwich  is  meant  a slab  of  oxidizer,  usually  ammonium  perchlo- 
rate (AP)  in  the  cited  investigations,  adjacent  to  a slab  of  polymeric  fuel. 
This  two-dimensional  configuration  is  an  important  experimental  configuration 
for  the  study  of  phenomena  taking  place  near  the  binder-oxidizer  interface 
during  a deflagration  process.  By  removing  the  difficulties  of  viewing  the 
interface  details  in  a real  three-dimensional  composite  solid  propellant, 
much  information  has  been  gathered  concerning  the  behavior  and  existence  of 
binder  melts,  catalytic  additives,  the  relative  importance  of  AP-binder  reac- 
tions in  driving  the  deflagration  rate,  fluid  mechanical  instabilities  during 
deflagration,  and  the  behavior  of  aluminum  when  placed  into  the  sandwich. 

In  contrast  to  the  many  experimental  studies  there  have  appeared  no 
ansl-ysea  of  the  sandwich  deflagration  process  with  AP  oxidizer  which  had  as  & 
goal  the  prediction  of  the  actual  shape  of  the  gas-condensed  phase  interface 
Bakhman  and  Librovich'^^  theoretically  investigated  a semi-infinite  slab  of 
oxidizer  adjacent  to  a semi-infinite  slab  of  fuel,  but  the  oxidizer  was  not 
assumed  to  be  capable  of  self-deflagration.  Nachbai investigated  a peri- 
odic, two-dimensional  array  of  oxidizer  and  fuel  slabs  as  a model  of  an  act- 
ual propellant.  The  goal  was  to  predict  a mean  deflagration  rate  without  pre- 
dicting an  actual  achieved  two-dimensional  surface  shape.  Again,  the  oxidizer 
was  assumed  incapable  of  self-deflagration. 

The  need  for  an  analytical  model  capable  of  prediction  of  the  sxu:face 
shape  for  an  AP  oxidized  sandwich  becomes  apparent  when  trying  to  reason  the 
behavior  which  is  experimentally  observed.  Even  in  this  two-dimensional  con- 
figuration the  problem  is  highly  complex  due  to  a)  the  appearance  of  binder 
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nelts,  b)  three  phase  heat  transfer  and  one  phase  (perhaps  two  phase)  mass 
transfer,  e)  two-dimensional  heat  and  mass  transfer,  d)  non-linearities  in 
the  governing  equations  due  to  at  least  ehetnieal  reaction  and  an  unknown  sur- 
faee  shape,  e)  multiple  ehemieal  reactions  and  condensed  phase  reactions, 
f)  a mathematically  elliptic  problem  in  an  infinite  region  and  g)  the  appear- 
ance of  an  eigenvalue  - the  deflagration  rate.  Under  these  difficulties  an 
attea^jt  to  reason  on  a qualitative  basis,  given  experimental  information,  is 
hampered  by  the  lack  of  some  con5)utational  framework  and  predictive  capa- 
bility. 

A previous  initial  step  toward  providing  a theory  of  the  sandwich  de- 
flq') 

flagration  process'  failed  to  yield  the  burn  rate  as  an  eigenvalue  of  the 
problem.  Purthermore,  an  analytical  error  has  been  discovered  in  that  treat- 
ment and  the  results  are  consequently  in  error.  The  current  treatment  cor- 
rects the  error  and  recovers  the  burn  rate  as  an  eigenvalue  of  the  problem 
by  a more  complex  solution  of  the  problem. 

Prom  a practical  standpoint  there  are  two  major  reasons  for  treating 
this  problem.  The  first  is  to  obtain  an  idea  of  the  distance  scales  in- 
volved in  the  problem.  Oiat  is,  how  far  do  the  influence  of  the  binder  and 


oxidizer  penetrate  into  each  other?  If  these  distance  scales  are  small 
enough,  this  interface  theory  may  be  xised  in  theory  of  an  actual  propellant. 
Secondly,  the  actual  burn  rate  is  desired.  Is  the  influence  of  the  binder 
on  the  self-deflagrating  AP  strong  or  weak?  Ihe  current  analysis  addresses 
those  issues. 


i 


37 


ANALYSIS 


Model  Construction  and  Assumptions 

Given  the  coraplicated  nature  of  the  problem,  a model  is  first  sought 
which  uses  available  ejtpcrimental  information  liberally  but  •irtiich  still  does 
not  overly  restrict  the  model  in  interpretation  of  experimental  results. 
Accordingly,  the  initital  model  uses  the  following  observations: 

a)  Far  from  the  binder-oxidizer  interface  the  AP  regresses  as  pure  AP. 
Furthermore,  for  binder  thicknesses  of  the  order  used  in  some  of  the  ex- 
perimental studies  (w  150  ^un)  there  is  little  effect  of  one  side  of  a sand- 
wich upon  the  other  side  even  when  dissimilar  materials  are  used.  IQiere- 
fore,  the  initial  model  development  is  concerned  with  a semi-infinite  slab 
of  AP  against  a semi-infinite  slab  of  binder. 

b)  A steady  state  is  achieved  experimentally  with  AP  oxidizer.  Conse- 
quencly,  time  dependence  is  assumed  absent. 

c)  For  uncatalyzed  sandwiches  the  e:q)erimental  results  show  very  little 
effect  of  the  binder-oxidizer  reactions  upon  the  surface  profile-  That  is, 
the  heat  feedback  from  the  binder-oxidizer  reactions  does  not  drive  the  over- 
all deflagration  rate;  the  AP  self- deflagration  is  responsible  for  the  over- 
all deflagration  rate  at  pressures  removed  from  the  low  pressure  deflagra- 
tion limit  and  below  2000  psia.  Hie  initial  model  is  therefore  constructed 
assuming  binder-oxidizer  reactions  to  have  negligible  rate.  This  does  not 
mean  that  a diffusion  flame  between  the  binder  and  oxidizer  is  absent,  but 
it  means  this  flame  is  not  close  enough  to  the  interface  to  play  a role  in 
the  deflagx'ation  behavior.  Furthermore,  the  effect  of  catalysis  is  not  treat- 
ed. 


_ 


The  initial  model  therefore  asks  the  question  of  the  surface  shape 
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attained  by  a semi-infinite  slab  of  AP  vhich  pyrolyzes  a f ^mi-infinite  slab 
of  binder.  Posed  in  this  manner  it  is  immediately  recognized  that  the  prob- 
lem has  neither  a unique  solution  nor  a soeady  solution  because  a)  the  final 
shape  would  depend  upon  the  geometry  of  ignition  and  b)  it  would  take  an  in- 
finite time  to  establish  a steady  profile  in  a semi-infinite  slab  of  inert 
binder.  The  ignition  problem  is  seen  by  imagining  two  cases  - one  in  which 
ignition  is  achieved  by  a line  heat  source  (say  an  ignition  wire)  and  a sec- 
ond in  which  ignition  is  achieved  uniformly  over  the  entire  AP  surface.  In 
the  first  case  the  AP  would  take  on  the  shape  of  ever-increasing  circular 
radii  from  the  ignition  point.  In  the  second  case  the  AP  would  deflagrate 
in  a planar  fashion  except  in  the  vicinity  of  the  binder.  However,  these 
comments  neglect  the  fact  that  the  inert  nature  of  the  binder  may  affect  the 
shape  attained  by  the  AP,  even  far  away  from  the  binder.  Surely,  at  long 
enough  time  the  profile  in  the  interface  vicinity  will  have  become  shape- 
invariant,  but  it  is  not  obvious  that  the  AP  will  be  horizontal  (assuming 
top-to-bottom  deflagration)  . In  fact,  this  is  to  be  determined  by  the  solu- 
tion. Concerning  the  binder,  it  appears  clear  that  after  a sufficiently  long 
time  the  processes  in  the  interface  vicinity  will  not  depend  upon  processes 
taking  place  at  a large  vertical  distance  from  the  interface.  Consequently, 
local  steadiness  may  be  presumed,  as,  in  fact,  e:q3erimentally  occurs. 

For  this  irdtisLl  model  the  absence  of  binder  melts  will  be  assumed.  The 
limits  of  validity  will  then  be  determined  by  con5>arison  of  the  model  and 

e^qierimental  results.  For  the  AP  deflagration  process  the  Guirao-Williams 

(ll) 

model  ' is  accepted  with  an  equilibrium  assuiqjtion  for  the  gas-solid  (or 
viscous  liquid)  AP  interface.  Some  minor  modifications  arc  introduced  into 
this  model  ['ir  conqjutational  convenience;  these  will  be  described  below.  Use 
of  this  model  will  restrict  the  sandwich  theory  vaLLidity  to  the  pressure 

L 


range  20-100  atm,  because  there  is  no  AP  theory  capable  of  an  explanation  of 
obcerved  jihcnoraena  abcjvc  100  aim  and  Ihc  low  pressure  deflagration  limit  of 
AP  occurs  near  20  atm. 

Other  usual  assumptions  are  made  to  simplify  the  analysis  which,  while 
they  lead  to  numerical  error  of  order  unity,  do  not  alter  significantly  the 
scaling  rules  developed  with  respect  to  other  variables.  These  assufl5>tions 
are:  a)  the  thermal  and  transport  processes  of  the  solid  AP  and  binder  are 

identical,  b)  ..he  thermal  and  transport  properties  of  all  gas  phase  species 
are  identical,  c)  the  Lewis  number  is  everywhere  unity  in  the  gas  phase, 
d)  the  deflagration  process  takes  place  at  constant  pressure,  e)  heat  con- 
duction and  mass  transfer  take  place  by  teoroerature  and  concentration  gra- 
dients, only,  respectively,  and  the  transport  coefficients  are  independent 
of  temperatxu'e  in  both  the  solid  and  gas  phases . A final  major  assumption 
is  that  on  any  vertical  line  parallel  to  the  binder- oxidizer  interface  the 
pv  product  (density  times  velocity)  is  that  as  determined  in  the  solid  phase 
and  all  lateral  velocities  are  zero  (strictly  true  in  the  solid  phase) . !Ehii 
is  in  the  spirit  of  the  Burke-Schumann  approximation  as  expounded  in  Refer- 
ence (2) . This  does  yield  error  in  convection  effects  \Q>on  heat  traxisfer, 
but  exact  treatment  of  the  problem  appears  too  complex  at  the  present  time. 

The  configuration  is  shown  in  Figure  1,  in  which  the  coordinate  system 
is  rendered  stationary  by  a translation  of  the  interface  in  the  y direction 
at  the  rate  r.  Under  the  stated  assumptions  the  equations  for  solution  and 
the  boundary  conditions  are: 

Gas  Phase  Species  Continuity 


* , * 
. ,'V”3 

* * 

T continuous,  VT  continuous  within  a phase 
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(U) 


(12) 


The  products  of  AP  gasification  are  assuned  NH^  and  HClOj^  which  are  assumed 

* 

identical  molecxiles  for  mass  transfer  computation,  k is  a rate  coef- 
ficient for  the  assumed  second  order  reaction;  as  written  in  Equation  (3) 

k is  not  a fundamental  preejqponential  constant  but  already  has  molecular 

*2 

constants  and  a factor  p absorbed  into  it.  The  factor  2 in  front  of 
q„w^  in  Equation  (4)  occurs  because  q^  will  be  quoted  per  unit  mass  of  AP 
rather  than  per  unit  mass  NI^.  The  equilibrium  interface  on  AP  is  spe- 
cified through  Equation  (lO)  . Uie  formulation,  as  far  as  the  AP  defla- 
gration process  is  concerned,  differs  from  that  of  Reference  (11)  in  the 
following  respects:  a)  no  dilution  of  the  and  HClOj^  is  assumed  at  tht 

solid-giiS  interface,  although  it  is  tacitly  accounted  for  by  the  choice 
* 

of  a number  for  qg.  calculations  are  singplified  by  taking  the  molecular 

weight  of  all  species  to  be  the  same.  The  constants  k and  b„  will  be  so 

J? 

chosen  to  recover  the  same  burn  rate  and  surface  tenrperat\ire  results  as 
in  Reference  (10) . 

Equation  (7)  is  the  pyrolysis  law  for  the  binder,  Equation  (8)  is 

the  energy  conservation  law  at  tht  solid-gas  interface,  auid  Equations 

(9)  and  (10)  are  the  interface  diffusion  laws.  Note  in  Equation  (8)  that 
* 

q undergoes  a discontinuity  at  the  binder-oxidizer  interface  and  Equa- 
tions  (9)  and  (10)  are  only  valid  on  the  AP  side  of  the  interface.  Shown 
in  Tabic  1 are  typical  values  used  in  this  work  for  the  various  parameters. 
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Table  1 

Kutacrlcal  Values  for  Various  Parameters 


Qyaiitity 


Value 


Reference 


> 

^0 

300  °K 

Assumed 

* 

o» 

1.95  gm/cm^ 

11 

=P* 

.3  cal/gm^K 

11 

< 

.3  cal/gm°K 

11 

173  oal/gm 

”to  yield  flf 
temperature 
1205°K  of 
Ref.  (11) 

* 

V 

-100  cal/gm 

11 

10“^^  cal/cm  sec°K 

* 

11 

9 X 10  cal/cm  sec^ 

il 

g 

15  kcal/mole 

11 

E* 

s 

30  kcal/mole 

11 

B 

8.7  - 17  kcal/mole 

14 

^B 

1 - 1;^  cra/sec 

14 

160  - 1004  cal/gm 

14 

“B 

Mathematical  Character  of  the  Problem 

Equations  (l)  - (5)  define  an  elliptic  problem  in  the  sense  thai  what 

happens  at  one  point  in  the  field  affects  every  other  point.  Further- 

* * * * 

more  an  eigenvalue  appears,  p v = p r . 

s 

If  the  binder  is  hard  to  decompose  and  it  Eissumes  a near’.y  vertical 
surface,  it  appears  obvious  that  the  picture  becomes  one  of  a (nearly) 
flat  plate  of  binder  ever  which  hot  AP  gases  are  flowing.  If  the  Reynolds 
number  based  on  distance  along  the  binder  were  large  enough  this  would 
revert  to  a parabolic  problem  because  3/3x  » d/dy  would  result.  However, 


exactly  at  the  binder-oxidizer  interface,  the  Reynolds  number  is  zero. 
Since  it  is  precisely  this  region  that  is  of  interest,  the  full  ellip- 
tic problem  must  be  solved.  In  order  to  gain  an  idea  of  magnitudes  in- 
volved here  the  equations  are  nondimcnsionalizcd  ^vith  respect  to  a dis- 
tance scale  and  temperature  T^.  Heats  of  gasification  are  made 

dimensionless  by  c*T*  and  activation  energies  by  R T^.  In  order  to  lo- 
cate the  condensed  phase  - gas  phase  interface  at  a constant  position 
the  y variable  is  replaced  by  T]  = y - yg(x)  . The  resulting  dimension- 
less equations  and  bcvindary  conditions  are 


ih]  - ° 

(Gas)  l(g]  = y/  e'V® 

(Solid)  l\s]  = (1  “ I)  g- 

Yp(x,»)  = 0 Yj.(-<»,y)  = Yj.(y) 

g(x,-=)  = 1 g(~“,y)  = g(y) 

g(oQ,y)  = 1 g(x,»)  = 1 


ag/dn(g^  - 


aYg/anl 


s+ 


= % 


AP 


-e  /g 

V ® 


s 


g continuous,  7g  continuous  within  a phaie 


(13) 


Since  y'  and  y*'  are  unknowns,  the  nonlinear  character  in  Equations (I3)  is 
s 

apparent.  Nonlinearities  also  arise  from  the  chemical  reaction  terms. 

Dio  dJmcnsion  a /r  is  known  to  be  the  "thickness"  of  the  thermal  wave 
s 

which  would  occur  in  a planar  regression.  It  is  the  reference  dimension 
here.  Die  parameter  5 in  Equations  (12)  is  nothing  more  than  the  ratio 
of  a characteristic  solid  phase  dimension  (ff  /r  ) to  the  characteristic 

S 

gas  phase  dimension  (or  /v  ) . If  there  were  no  modification  due  to  the 
reaction  rate  term,  the  gas  phase  distance  over  which  significant  heat 
transfer  would  occur  would  be  of  the  order  of  a/v  . Using  the  para- 
meters of  Table  1,  ? = 9*0,  showing  that  the  gas  and  solid  phase  charac- 
teristic scales  arc  quite  different.  Furthermore,  constructing  the 
Reynolds  number  based  upon  y,  it  is  found  that  Re  = %y,  so  that  when 

V 

y is  of  the  order  of  1/5  a transition  is  taking  place  between  "low" 
and  "high"  Reynolds  numbers.  If  important  field  quantity  variations 
are  taking  place  only  over  a gas  phase  distance  of  the  order  of  1/5, 
the  problem  must  be  treated  as  elliptic  with  no  simplifications  possible 
through  a boundaiy  layer  assumption.  However,  a rather  simple  method 
of  solution  will  emerge  if  at  some  point  in  the  flow  field  the  boundary 
layer  approximation  may  be  involved.  The  use  of  this  approximation  will 
be  illustrated  below. 

Solution  by  an  Integral  Technique 
Pure  AF  Deflagration 

Far  from  the  binder  the  AP  must  undergo  a planar  deflagration  but 


the  angle  9 is  unknown;  it  is  the  eigenvalue  of  the  problem.  All 


U5 

x-derivativcs  must  vanish  and  Equations  (I3)  hecome  nonlinear  ordinary 
differential  equations  in  T\  along  vrith  the  appropriate  boundary  condi- 
tions. 'ihc  equation  for  the  solid  phase  heat  transfer  may  b>»  solved 
exactly.  'Ihe  solution  is 

2 

i “ 1 » (ig  - 1)  (1^) 

An  overall  energy  balance  yields  the  adiabatic  flame  tei^ierature 

Sf  = «s 

A first  integral  of  the  gas  phase  heat  and  mass  transfer  equations, 
subject  to  the  boundary  conditions,  is 

a(Tt)  = g(lO  + 2qpyp,('fl)  « g^  = (I6) 

s 

Ihe  solut'  on  is  cojjpleted  by  assuming  a functional  form  for  g(T))  as 


g(lD  - Kg  = (Sf 

- ig)  ? 
c 

11  < c 

g(Tl)  = gf 

11  i c 

(17) 

>rtiere  c is  the  11  position  where  y„(T0  vanishes,  i.e.  reaction  is  cotaplete. 

X 

Now  integrating  the  gas  phase  energy  equation  between  11  = 0 and  T1  * c 
using  the  aasiuiied  form  for  g(ll)  from  Equation  (16)  , using  the  solid-gas 
interface  boundary  conditions  and  the  equilibrium  condition  from  Equations 
(12) , the  solution  for  the  AP  deflagration  is  completed  as 


K/S  ; (l  - ‘ d(  2 ) 


Yp  (/+  ?)  = 5/2 


y = e 
F^  F 


■®S  /®8 


An  extremely  in5>ortant  point  to  note  from  Eqs.  (l8)  is  that  they  ing>ly 
that  the  regression  rate  normal  to  the  gas-solid  interface  is  independent 
of  9,  the  inclination  to  the  horizontal.  For  note,  Q(c)  is  an  invariant. 

For  a fixed  g , is  fixed  and  c “ z^.  Therefore,  Q “ kY?  c « kz^.  But 
k “ 1/t  so  that  r « z and  r = r cos  0 = r/z  is  an  invariant.  Conse- 
quently,  for  any  sandwich  vertical  regression  rate,  r,  which  is  the  eigen- 
value of  the  problem,  the  AP  far  from  the  binder  will  regress  normal 
to  its  surface  at  a fixed,  unique  value.  The  pertinent  values  for  AP 
combustion  are  consequently  only  presented  for  0=0.  Tto  numerically  com- 
plete the  solution  the  following  procedure  is  used,  a)  as  a function 
of  pressure  is  taken  from  Reference  (11) ; b)  Equation  (l4)  yields  g^ 

(\diich  is  actually  constant  here  because  = 1 and  q^  and  q^  are  assumed 
independent  of  pressure) ; c)  Equation  (16)  determines  ; d)  Equations  (1?) 
detemune  o,  bp  and  k. 

Prom  the  nondimens ionalizati on  procedure  it  may  be  recalled  that  k « 

2 -2 

p /r  . Therefore,  if  the  rate,  r,  is  known  at  one  pressure,  it  is  known 
as  a function  of  prcssxire.  Shown  in  Table  2 are  conq)lete  calculations 
for  two  35Jt^le  pressures. 

In  Table  2 it  will  be  noted  that  b„  is  not  quite  constant.  This  is 

£ 

due  to  the  use  of  a slightly  higher  E than  in  Reference  (10)  . Furthermore, 

®F 
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k is  not  quite  proportional  to  (p/r)  . “Riis  is  due  to  the  fact  that  r 
in  Table  2 is  the  ejfperiracntal  value  and  it  is  known  that  between  the  two 
values  of  pi'essurc  of  Table  2 the  theory  of  Reference  (11)  slightly  over- 
cstiiaal-os  the  prcssxu’e  sensitivity  of  the  deflagration  rate.  These  de- 
tails are  not  considered  important  for  the  current  theory  because  a pre- 
cise model  for  AP  deflagration  is  not  sought;  only  the  deviations  from  a 
planar  regression,  due  to  the  binder  presence,  are  required  and  the  above 
theory  appears  adequate  to  serve  as  a baseline  for  perturbations  due  to 
the  binder. 


Table  2 


One  Dimensional  AP  Deflagration  Results 


p 

r 

(atm) 

(cm/sec) 

(|im) 

54.4 

.735 

20.9 

100 .0 

1.000 

15-3 

®s  ®f 

(°K) 

880  2.93  .285  4.022 
911  3-04  .258  4.022 


(atm) 

.1471  4.38  X 10® 

.1183  4.06  X 10® 


k 

2.62  X 10® 
3.39  X 10® 


Perturbed  Solution 

It  is  noted  that  there  is  only  a very  weak  variation  of  the  dimensionless 

parameters  with  pressure  in  Table  2 (because  g variations  are  weak  with 

s 

pressure)  and  therefore  Eqs.  (I3)  are  nearly  pressiire  invariant.  Conse- 
quently, all  further  work  will  be  carried  out  for  the  conditions  that 

g„  = 2.93  corresponding  to  p = 54.4  atm.  = 800  psia. 
s 

Upon  extensive  investigation  of  the  partial  differential  equations 
for  small  deviations  from  the  one-dimensional  regression  it  was  determined 
that  a)  the  deviations  from  the  planar  case  in  the  gas  phase  could  be  ex- 
pected to  be  simple  deviations  from  the  planar  solution,  but  b)  the  solid 
phase  deviations  roay  be  complex.  By  "sickle"  it  is  meant  that  the  de- 
viation is  not  oscillatory.  Thus,  if  g(x,y)  = g(x)  +6  (x,y)  where  G (x,y) 
is  the  deviation  from  the  pure  AP  case, G (x,y)  may  be  expected  to  have 


monatonic  behavior  in  7 between  the  two  end  values  G [x,yj.(x)  ] . Ihere- 
forc  it  vas  decided  to  utte.rpt  an  integral  solution  where 


~ [}  ' 


Tl  < c 


Y«Y  +(Y  -Y')'^ 

B ' B B c 

s c s 


Yp-0 


Ti  < c 


Tl  > c 


g - gg(x)  = [gi(x)  > gg(x)  ] ^ Tl  < c (19) 

are  guessed  forms  of  the  solution  for  the  gas  phase. 

If  the  Eqs.  (19)  are  placed  in  Eqs.  (13)  and  integrated  from 
Tl  = 0 to  c(x)  the  result  is  three  nonlinear  ordinary  differential  equa- 
tions in  tlie  unknowns  Yj,  , Y^  , Y^  , c,  g^,  y'  and  the  Tl  deriva- 

■'^s  s “c 

tives  of  g and  Y,.  at  Tl  = c.  These  derivatives  appear  because  no  func- 

D 

tional  form  is  specified  for  Tl  > c.  It  is  anticipated  that  these  de- 
rivatives will  be  very  small  beyond  the  reaction  region  euid  one  pos- 
sible assumption  is  that  they  are  zero.  An  alternate  assumption  has 
been  employed  and  is  dlicussed  below.  The  diffusion  boundary  conditions 
and  the  equilibrium  interface  boundary  condition  provide  three  more  rela- 
tions for  the  nine  unknowns.  At  1^  = c which  is  of  the  order  of  magni- 
tude of  l/§,  the  ReynoT-Ts  number  is  making  a transition  from  low  to  high 
values,  and  it  appears  at  this  point  reasonable  to  assume  that  — » 3/3y  = 

A 

9TI,  Making  this  assuu5)tion  in  Eqs.  (13)  there  results 

BTl  /_  ~ > 2 /.n  ~ I 1 2 

‘T|  = c ^ox  Tls=c  ^dx 


B N 1 ^ > 1 

“5ti'  /Tl  = c’  ? ^.2  y Tl  = c"  ? 


^ = dx^ 


(20) 


Sqs.  (20)  eliminate  two  of  the  above  unknowns  and,  one  further  relation 
is  needed.  For  algebraic  siagilicity  another  difierential  equation  is 
obtained  by  taking  an  T|  - moment  of  the  equation 

L[a]  =0  a = g + 2qpY^,j 

deriveable  from  Eqs.  (I3) > and  integrating  from  o to  c. 

In  order  for  a numerical  integration  to  proceed  smoothly  to  x * 
an  asymptotic  solution  is  desirable.  Consequently,  a solution  was  first 
sought  which  is  a small  perturbation  about  the  planar  AP  deflagration 
state.  Letting 


®1  ■ ®1 


c = c + C(x) 


Vg  - l(!t)  . 


Y = Y (x) 
s s 


\ = Vb  W 
c c 

substituting  these  forms  into  the  nonlinear  ordinary  differential  equations 
and  boundary  conditions,  making  use  of  the  AP  solution  properties,  and 
neglecting  prod’’cts  and  squares  of  perturbation  quantities,  there  results 
the  foil  owing  linearized  set  of  equations  and  boundary  conditions 

Ky,  +yV;  = Q G+CL  Y + C + Q G (2l) 
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Gj^(-“)  = (|(-“)  " \ (••“)  * Yg  (-“)  = Vg  (-“)  = Y"(““)  = C(-»)  = 0 (27) 


Here  Qo.  Q Q and  Qv  partial  derivatives  of  the  reaction  rate 

integral  which  may  he  ® numerically  evaluated.  Note  that  Eqs . (23) 
and  (24)  are  decoupled  from  the  rest  of  the  system.  Because  Eqs.  (21)- 
(27)  are  linear  and  homogeneous,  they  possess  solutions  like,  say,  => 
AG  e”^.  Since  2Yp  + Yg  + Y^  = 1,  the  binder  equations  must  also  have  so- 
lutions  Yg  - Ay  where  the  m is  the  same.  So  there  are  two  sets 

of  equations  vrtiicfi  may  be  investigated.  Since  the  equations  are  homo- 
geneous, the  determineint  of  the  coefficient  matrix  must  equal  zero  for 
a solution  to  exist.  From  Eqs.  (21),  (22),  and  (25)  m = m^(0)  is  there- 
fore developed.  From  Eqs.  (24)  and  (25)  m = 1112(0)  is  developed  and  the 
results  are  shown  on  Fig.  2.  Only  real  values  of  m were  numerically  in- 
vestigated because  the  only  unique  solutions  which  will  properly  attach 
to  the  binder  are  for  real  m,  as  will  be  seen  later.  It  is  seen  that 
there  is  a unique  9 for 

m = 6.4  0=2° 


That  is,  an  eigenvalue  exists.  The  AP  assumes  a nearly  horizontal  sur- 
face fur  from  the  binder  so  that  the  burn  rate  is  only  slightly  higher 
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and  the  physical  properties  of  thebinder.  A question  now  to  be  asked  is  whe- 
ther or  not  this  linear  solution  may  be  used  instead  of  a full  nonlinear 
solution  for  binders  of  practical  interest,  ihe  usual  physical  ex- 
pectation is  that  a linearized  theory  will  be  reasonable  if  the  pertur- 
bation quantities  remain  within  some  prescribed  fraction  of  the  cor- 
responding baseline  quantity.  Vi'iwing  Fig.  3 and  focusing  on  g^,  for 
example,  if  an  arbitrary  limit  of  validity  is  set  that  (g^  - g2^)/g^ 

< 10^,  it  appears  that  the  solution  might  be  expected  to  be  valid  for 
X <0.6.  Adopting  this  criterion,  a procedure  to  match  to  the  binder 
is  required  for  the  linearized  solution. 

Location  of  the  Binder 

It  will  be  noticed  that  the  above  solution  is  independent  of  the  so- 
lid phase  solution.  In  fact  the  heat  transfer  condition  of  Eqs.  (13) 
and  the  surface  temperature  of  the  gas  phase  solution  form  boundary 
conditions  for  the  solid  phase  heat  trsinsfer.  Since  the  attenqjt  here 
is  to  work  with  the  asymptotic  solution  to  a linear  set  of  equations , 
the  consistent  linearized  solid  phase  equation  from  Eqs.  (I3)  is  for 
S - iCTl)  + G(x,Tt) 

ax2  = ^ all 

where  the  sic^jlifVing  assurnption  of  y^=0  has  been  made  in  accordance  with 

s 

the  above  gas  phase  solution.  Letting  G = g + Y t§,  Eq.  (2?)  becomes 

oil 


which  has  the  solution 


0 


g = e^^  e*”^  [g  cos  m_  11  + G,  sin  m_  T|] 
o g 1 S 


(29) 


This  'lolut.lcn  has  a feature  that 


S=  = V“  * ^ II 

which  determines  from  the  gas  pheise  solution.  A most  important  feature 
of  the  solution  is  the  creation  of  an  oscillatory  ripple  of  very  short 
wavelength  (of  order  l/ra  w l/m)  superimposed  on  the  monotonicaJLly  vary- 
ing  temperature  of  the  undisturbed  solution.  What  happens,  therefore, 
is  that  in  both  the  solid  phase  and  gsis  phase  the  distance  scales  over 
which  rapid  transitions  are  made  are  now  of  the  same  order  of  magnitude 
(of  order  l/g)  . In  a sense,  the  short  distance  sceJ-e  of  the  gas  phase  is 
imposed  upon  the  solid  phase. 

TO  complete  the  solid  phase  solution  Eq.  (28)  is  substituted  into  the 
interface  heat  transfer  condition  of  Eqs.  (I3)  and  is  determined.  Now  at 
any  x position  in  the  eigensolution  to  the  AP  problem  the  heat  transfer  vec- 
tor in  the  gas  and  solid  phases  is  known.  Since  this  must  be  a contin- 
uous quantity  and  the  temperatxire  is  a continuous  quantity,  but  q under- 

s 

goes  a discontinuity,  there  must  be  a surface  slope  discontinuity  at  the 
binder.  In  the  interface  energy  conservation  relation  of  Equations  (I3) 
the  solid  and  gas  phase  heat  transfer  vectors  may  be  coniputed  fi*om  the  AP 
solution  and  this  equation  becomes  a relation  for  the  binder  heat  of  gasi- 
fication as  a function  of  its  surface  slope.  Ihe  result  is 


gi  - 


- (g^  - 1)  (1  + Y)  - ^ - m gJ  (30) 


The  pyrolysis  condition  of  Equations  (I3)  gives  an  additional  relation  be- 
tween the  surface  slope  and  the  binder  properties 
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1 ^ s 

- = b,e  ® 


At  any  x position,  then,  Equations  (29  and  30)  together  with  the  A? 
eigensoiution  define  an  allowable  binder  attachment  and  a functional 
equation 


(31) 


“ ^B 


4 = (%■%-,’=') 
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For  the  case  of  p = 54.4  atn  these  results  are  shown  in  Figures  4 and  5 for 

(14) 

two  values  of  e v^iiich  correspond  closely  to  HTPB  and  CH’B  binders. 

(14) 

Also,  knowing  q and  b_  for  HIPB  and  CITB  ‘ the  actual  point  at  which 

S D 

1 * 

these  two  binders  would  attach  is  shown  on  Figures  4 and  5.  Ohere  are 
several  points  worthy  of  note.  First,  for  these  binders,  the  x posin-ion 
of  con^jatability  with  the  AP  solution  occurs  where  very  little  change  from 
a flat  AP  surface  has  taken  place.  Consequently,  the  linear  AP  eigensolu- 
tioii  can  be  used  with  confidence  as  a good  approximation  to  the  solution 
of  the  nonlinear  problem.  Polyurethane,  shown  on  Fig.  4 would  also  at- 
tach to  a nearly  flat  AP  surface.  A check  of  the  fluorocarbon  binder 
data  of  Ref.  (l4)  also  shows  an  attachment  at  negative  x.  The  second 
imp(^rtant  point  is  that  the  binder  slope  is  very  nearly  vertical  at  the 
attachment  point.  This  is  the  primary  information  desired,  in  addition 
to  the  AP  surface  profile,  so  no  attempt  is  made  to  continue  the  solution 
to  tho  right  of  the  attachment  point  to  find  the  binder  profile.  The 
third  major  point,  refer! ng  back  to  Fig.  3»  is  that  the  solid  phase  heat 

■|1u»  activation  energies  for  HTPB  and  CTPB  do  not  exactly  fit  the 
numbers  of  Figs.  4 and  5»  What  has  been  done  is  to  compute  a b„  to  fit 
the  known  binder  pyrolysis  data  at  a teraperature  g , assuming  tne  acti- 
vation energies  of  the  figures.  ® 
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trausfoi-  vof^tor  points  from  the  binder  toward  the  AP.  This  indicates  that 
at  least  for  a short  distance  above  the  interface  the  binder  temperature 
must  bo  increeising,  probably  due  to  the  fact  that  the  highest  temperature 
AP  gases  are  at  positive  t). 


DISCUSSION  OF  PESUI'TS 


A solution  has  been  obtsdned  for  the  shape  of  the  deflagrating  AP 
surface  when  it  is  adjacent  to  an  inert  dry  binder.  Ihe  eigenvalue,  the 
vertical  regression  rate,  has  been  found  to  be  independent  of  binder  type 
8Uid  to  deviate  only  very  slightly  from  the  pure  AP  burn  rate,  as  has  been 
e3q>erimentally  observed.  Except  for  selected  binders  there  would  be  very 
little  visible  effect  of  the  binder  upon  the  surface  shape  and  the  result 
is  virtually  independent  of  pressure.  Ihe  distance  scale  over  which  a 
visible  transition  would  take  place  from  planar  AP  to  the  binder  is  of  the 
order  of  microns.  The  current  theory  assumes  a dry  binder;  it  is  known, 
however,  that  binder  melt  flows  exist  for  all  binders  tested  heretofore 
in  the  sandwich  configuration  and  that  these  melt  flows  run  several  hun- 
dred microns  onto  the  AP  surface.  Furthermore,  surface  roughness  demen- 
sions  of  the  order  of  micronsdevelops  during  deflagrations.  Consequently, 

• . I 

few  of  the  predicted  phenomena  are  capable  of  being  observed.  A theory 
including  the  effects  of  melt  flows  is  necessary. 

The  theory  predicts,  however,  that  if  the  melts  do  not  occur  there 
should  bo  a sharp  discontinuity  in  slope  at  the  binder-oxidizer  interface. 
This  has  recently  been  seen  for  catalyzed  sandwiches for  which  the  melt 
extent  is  markedly  reduced  (for  unknown  reasons) . Although  the  current 
theory'  is  not  directly  applicabl3  to  catalyzed  situations  the  interface 
conditions  responsible  for  the  slope  discontinuity  arc  applicable.  It 
appcfu‘6  that  melt  flows  dominate  the  development  of  the 


surface  shape  if  melcs  occur. 

Iho  current  theo  t shows  that  the  heat  flow  vector  near  the  gas- 
soli'l  interface  and  Al-binder  solid  interface  is  from  the  binder  toward 
-no  AJ’  in  the  solid  i>hase.  Consequently,  there  must  exist  a weak  (hot) 
portion  of  the  binder  slightly  up  from  the  four-corner  interface  and  into 
the  binder.  This  might  account  for  the  appearance  cf  "notches"  in  the 
binder  sometimes  seen  in  quenched  saciples . The  violence  of  the 
quenching  process  may  eject  the  part  of  the  binder  which  is  weaker  than 
surrounding  parts. 

The  current  analysis  shows  a very  weak  dependence  of  surface  shape 
vroon  pressure.  This  independence  has  been  observed  experimentally,  but 
the  comparison  between  theory  an]  experiment  cannot  be  made  precisely  be- 
cause binder  melts  have  occurred  in  sJJ.  the  ekperimenta. 

If  the  AP  particle  size  is  sufficiently  large  the  current  euiadysis 
may  form  the  basis  for  a deflagration  theory  of  a heterogeneous  propellant, 
if  the  binder  were  dry.  The  largest  natural  dimension  which  occurs  in 
this  theory  is  the  thermal  wave  depth  which  is  of  the  order  o*"  co  pm  at 
800  psia.  Consequently,  for  AP  particle  sizes  larger  than  this  dimension 
there  maybe  some  merit  in  applying  this  technique  to  a real  propel-Lant, 

An  extension  to  the  case  of  a finite  binder  thickness  would  be  required, 
however,  unless  AP  particle  sizes  substantially  in  excess  of  the  thermal 
wave  depth  were  considered  because  the  typical  binder  widths  would  be  less 
than  the  typical  thermal  depth  in  the  binder.  Such  an  extension  is  not 
deemed  difficult  to  attain. 
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GAS  PHASE 


Figure  1.  Sandwich  Sche;-natic  and  the  Coordinate  System 


Figure  3.  Surface,  Plane  Standoff,  Surface  TJIL  Mass  Fraction,  Flame 
Temperature  and  Solid  Phase  Heat  Trrbsfer  Vector  Profiles 


